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ABSTRACT 

The observation of GRB 0803 19B, with an isotropic energy E iso = 1.32 x 10 54 erg, and GRB 050904, with 
Eiso = 1.04 x 10 54 erg, offers the possibility of studying the spectral properties of the prompt radiation of two 
of the most energetic Gamma Ray Bursts (GRBs). This allows us to probe the validity of the fireshell model 
for GRBs beyond 10 54 erg, well outside the energy range where it has been successfully tested up to now 
(10 49 -10 53 erg). We find that in the low energy region, the prompt emission spectra observed by Swift BAT 
reveals more power than theoretically predicted. The opportunities offered by these observations to improve 
the fireshell model are outlined in this paper. One of the distinguishing features of the fireshell model is 
that it relates the observed GRB spectra to the spectrum in the comoving frame of the fireshell. Originally, 
a fully radiative condition and a comoving thermal spectrum were adopted. An additional power-law in the 
comoving thermal spectrum is required due to the discrepancy of the theoretical and observed light curves 
and spectra in the fireshell model for GRBs 0803 19B and 050904. A new phenomenological parameter a is 
correspondingly introduced in the model. We perform numerical simulations of the prompt emission in the 
Swift BAT bandpass by assuming different values of a within the fireshell model. We compare them with the 
GRB 0803 19B and GRB 050904 observed time-resolved spectra, as well as with their time-integrated spectra 
and light curves. Although GRB 0803 19B and GRB 050904 are at very different redshifts (z=0.937 and 
z=6.29 respectively), a value of a = —1.8 leads for both of them to a good agreement between the numerical 
simulations and the observed BAT light curves, time-resolved and time-integrated spectra. Such a modified 
spectrum is also consistent with the observations of previously analyzed less energetic GRBs and reasons for 
this additional agreement are given. Perspectives for future low energy missions are outlined. 

Subject headings: Gamma-ray burst: general — Gamma-ray burst: individual: GRB 0803 19B — Gamma-ray 
burst: individual: GRB 050904 — ISM: structure — Black hole physics 



1. INTRODUCTION 

Out of the hundreds of Gamma Ray Bursts (GRBs) so far 
observed with known redshifts, there are approximately ten 
GRBs with an isotropic energy E iso > 10 54 erg: GRB 990123, 
GRB 990506, GRB 000131, GRB 050820A, GRB 050904, 
GRB 0803 19B, GRB 080607, GRB 0807 21, GRB 0809 16C 
GRB 090323, GRB 09 0926A (see e.g. lAmati et alj 120081 
20091: iKann etalJ 120101) . GRB 090902B (see jAbdo et alJ 
2009) and GRB 110918A (Eederiks & Pal'ShinlBOl ll) . We 



will analyse two of these sources to probe the fireshell model, 
which has been successfully applied to GRBs with Ei S0 up 
to 10 53 erg. The two candidates are GRB 0803 19B and 
GRB 050904, having an isotropic 7-ray energy release re- 
spectively of E iso = 1 32 x 10 54 erg (20 keV - 7 MeV, see 
iGolenetskii et al J2008I) and E,,„ = 1 .04^5 x 10 54 erg (15 ke v 
- 5 MeV. see lSugita et al.ll2009l) . 

Much of the progress made in observing GRBs in recent 
years has been due to the coordinated efforts of a large num- 
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ber of satellites including Ko nus- WIND ( Apt ekar et al.l 19951). 
Swift dGehrels et al.l 120041) Suzaku-WAM dYamaoka et al] 
2009), AGILE dTavani et"al]l2009l) and Fermi (lAtwood et all 
2009; Meegan et al. 2009). These satellites allow an overall 
energy coverage from 0.2 keV to 300 GeV. 

GRB 0803 19B and GRB 050904 have been triggered by 
Swift. What is relevant here is that, in both these sources with 
unusually high Eiso, the observed peak energy E°^ k occurs 
well above the Swift BAT bandpass (see Tab. [3] in Sec. |6]): 
for GRB 050904 we have E°£ s ak = 314!™ keV iSugita etail 
( 120091) ; for GRB 0803 19B we have E° b e s ak = 675 ±22 keV. 
This, in turn, through the Swift observations, for the first time 
allows the exploration of GRB spectra at E < 0.1E peak (see 
Fig. [TJ. This is at variance with the previous observations 
of lower energetics GRBs where E peak falls within the in- 
strumental bandpass (see Fig. [TJ. The observation of GRB 
0803 19B has occurred prior to the launch of Fermi and after 
the one of AGILE, but AGILE has been occulted by Earth dur- 
ing the burst detection. Th e high energy ph otons have been 
detected by Konus-WIND dRacusin et all 120081) . The obser- 
vation of GRB 050904 occurred prior to the launch of both 
AGILE and Fermi. The high energ y photons were de tected 
by Suzaku-WAM and Konus-WIND (Sugit a et al.ll2009l) . 

GR B 0803 19B was dis covered by BAT on March 19, 
2008 dRacusin et all l2008h . It has a redshift z = 0.937 
dVreeswiik et al.ll2008l)~ and is characterized by an extraordi- 
narily bright optical emission accompanying its 7-ray emis- 
sion, that makes it the brightest optical burst ever observed: 
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with a peak visual magnitude of 5.3, it could have been 
seen with the naked eye by an observer in a dark location 
dRacusinet all 12008). It is also one of the brightest GRBs 
in 7 and X rays. 

This extremely bright optical flash has promoted an intense 
theoretical analysis on the role of synchrotron self-Compton 
(SSC, see Sec. |2]l radiation and this scen ario has been investi- 
gat ed by several authors, includ ing e.g. iRacusin et al.l (120081) 
and lKumar & Panaitescul ( 120081) . 

There has been also interest in examining the possible ex- 
istence of beaming in order to reduce the energetics of this 
source. IRacusin et al.l (12008) proposed a two-component jet 
model: a narrow jet, with a half-opening angle 0j et ~ 0.2°, 
dominating the early emission, surrounded by a wider jet, 
with 9 j et ~ 4°, dominating the emission at intermediate and 
late times. This two-jet model would reduce the total energy 
budget to ~ 4 x 10 50 erg. However, the narrow jet should 
produce a jet break ~ 1 hr po st-burst, which has not been 
observed. IRacusin et al.l ( 2008) suggested that this could be 
explained if the optical flux of the narrow jet is fainter than 
that of the wider jet. The wider jet should give rise to a break 
at ~ 10 6 s. iTanvir etal] ([2010) reported the observation of a 
jet break at ~ 1 1 days with an "approximately achromatic be- 
haviour", recalling that, in the Swift era, the expected "achro- 
matic behaviour of breaks in X-ray and optical light curves 
has been rarely seen (ICurran et al.ll2008l)". An alternative sin- 
gle jet model has been considered by iKumar & Panaitescul 
(2008): taking into account the lack of an optical jet break 
during the first ten days of emission, they find a lower limit 
6 j e t > 2°. The energy of the GRB could in this case be 
> jo 52 - 3 erg (twice larger for a double-sided jet). 

GRB 050904 was dis covered by BAT on September 4, 2005 
(Cummings et al. 20QH). It is one of the fa rthest GRBs ever 
observed, with z = 6.29 dKawai et alJl2005l) . Also this burst 
is characterized by an intense optical emission: a bright opti- 
cal flare was in fact detected by TAROT near the end of the 
prompt phase, in temporal coincide nce with an X-ray flare 
(Boe r et al.ll2006tlGendre et al.ll2007l) . 

Tagliafer ri et al.1 (120051) found a "steepening" in the /-band 
light curve at 2.6 ± 1.0 days and proposed that "it may be due 
to a jet break". In this case a jet opening a ngle Q\ Kt ~ 3° coul d 
be inferred. This analysis was refined by iKann et ail (12007b . 
who put the steepening at 2.63 ± 0.37 days, corresponding 
to 9 jet ~ 3.30°. In this case the beaming-corrected energy is 
1.73 x 10 51 erg dSueita et al. 1120091) . 

In this paper we limit the analysis of these two sources to 
the 7-ray emission of the prompt phase, which is energetically 
predominant with respect to the opt ical emission: t he optical 
isotropic energy is ~ 0.1% of Ej SO (Bartolini et al. 2009) for 
GRB 0803 19B and is eve n less for GRB 050904 dBoer et alJ 
120061 lGendre"eraIll2007h . 

Also in view of the absence of achromatic breaks required 
by the jetted emission model, for GRB 0803 19B and GRB 
050904 we assume spherical symmetry, which is one of the 
main features of the fireshell model. 

We recall that a satisfactory agreement between the the- 
oretical predictions of the fireshell model and the observed 
light curves and spectra of GRBs with Ej m up to 10 53 erg 
has been previously obtained (see e .g. Berna rdini et al.ll2005t 
Ruffinietai] 120061: iDainotti etal] 120071: iBernardini et all 
20071: iCaito et al.ll2009L I2010L Ide Barros et al.11201 11) . What 
is really new in the analysis of these two sources is: 1) the 
very large value of E° b e s ak , well above the Swift BAT bandpass 
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Figure 1. Theoretically simulated spectra obtained by assuming two differ- 
ent values of £, JO : 10 54 erg (in red) and 2.5 X 10 53 erg (in black). Solid lines 
refer to the case a = —1.8. Dotted lines refer to the ca se a = 0.0, correspond- 
ing to the pure thermal comoving spectrum (see Sec. 13.3.21 . The box marks 
the energy band covered by Swift BAT (15 keV - 150 keV). The behaviour of 
the spectra ar ound th e peak energy is not significantly affected by the value 
of a (see Sec. 13.3.21 . For the higher energetic sources, the low energy com- 
ponent of the spectra falls inside the energy band covered by BAT; this is not 
the case for the lower energetic sources, for which only the region around the 
peak can be investigated. 

(see Fig. [T]and Tab. [5); 2) the fact that for both sources the 
Swift BAT bandpass covers the low energy part of the spec- 
trum, that could then be investigated; 3) in the case of lower 
energetic sources, the energy region near the peak has been 
observed by BAT, while the low energy spectral component 
has been missed (see Fig.Q]) and will possibly be the object of 
future specific space missions. 

In both GRB 0803 19B and GRB 050904, Swift BAT data 
show more power in the low energy region than the theo- 
retical predictions of the fireshell model (see Sec. |4). We 
then introduce a new phenomenological spectrum in the co- 
moving frame o f the f ireshell: the "modified" thermal spec- 
trum (see Sec. I3.3.21 i. characterized by a phenomenologi- 
cal parameter a. The new comoving spectrum: 1) allows 
the correct reproduction of the observed prompt 7-ray emis- 
sion light curves and spectra of both GRBs (see Sees. [4] and 
[5]); 2) clearly, as shown by Fig. Q] does not modify signifi- 
cantly the considerations on the less energetic GRBs previ- 
ously discussed within the fireshell model, since in that case 
only the region encompassing the peak is inside the instru- 
ment bandpass and is not significantly affected by the value 
of a; 3) predicts a possible broader emission in the low en- 
ergy spectral component of the lower energetics GRBs which 
will possibly be tested by future missions below 10 keV such 
as, e.g., LOFT (Feroci & The LOFT Consortium! 1201 II) and 
MIRAX (lAmati etalJ 120111) . These observations will de- 
termine the possible general validity of the comoving spec- 
trum introduced here. It will also be important to verify if 
an analogous phenomenological correction in the Wien part 
of the comoving blackbody spectrum will allow the inter- 
pretation of the high energy data (above hundreds of MeV) 
later observed in equally energetic GRBs by AGILE and 
Fermi. In this respect, it must be noted that an extra high- 
energy power-law component has been actually invoked in 
the current literature to explain the emissi on above hundred s 
MeV observed in e.g. GRB 090902B dAbdo et alJ 2009), 
GRB 090510 (AckermanniLiDISOKl) and GRB 090926A 
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dAcker mann et al.l201 lb . However, the phenomenological ex- 
tra power-law component characterizing the low energy tail of 
the new comoving spectrum introduced here is not necessarily 
related to the low energy extrapolation of such a high-energy 
component. In fact it is also present in sources which show 
no evidence of h igh energy emissio n, like e.g. GRB 09061 8 
dlzzo et al.ll20T2b and GRB 101023 dPenacchi oni et al. 20p). 

The work is organized as follows. In Sec. [2] we describe 
the phenomenological and theoretical interpretation of GRB 
prompt emission. In Sec. [3] we briefly outline the fireshell 
model, we introduce the "modified" thermal spectrum and ex- 
plain how we perform the numerical simulations to be com- 
pared with the observational data. In sees. |4]and[5]we present, 
as specific examples, the analysis of GRB 0803 19B and GRB 
050904 observations of the prompt 7-ray emission. In Sec. [6] 
we discuss the possible explanation for the need to introduce 
the "modified" thermal spectrum and in Sec.|7]we present our 
conclusions. 

2. THE SPECTRAL ENERGY DISTRIBUTION OF GRB PROMPT 
EMISSION 

The study of the GRB prompt emission is often per- 
formed by means of phenomenological approaches consist- 
ing, for example, in modeling the burst spectra with empirical 
functions whose parameters are fixed by the \ 2 minimiza- 
tion procedure of the spectral da t a (see, e.g., jPreece et al.l 
1 19981 120001 fGhirlanda et al.M2002t iKaiieko et al.ll2006L 120081 
iGuidorzi et al.112.01 lb . The explanation of the results obtained 
from such phenomenological approaches are a main goal of 
the theoretical models of GRBs. 

2.1. Phenomenological approaches to GRB prompt emission 

The photon number spectrum time-integrated over the GRB 
prompt emission duration is typically best fit by two power- 
laws jomed_smorthly_at a given break energy (the "Band func- 
tion", iBand et aTll 19931) . whose low energy and high energy 
photon indic es, a and /3, have median values of — 1 and -2.3 
respectively dPreece et al. 2000; Kaneko et al. 20061). 

Another example of phenomenological analysis is the cor- 
relation between the rest-frame value of the peak energy E' peak 
and the isotropic equivalent radiated 7-ray energy 2? !SO , which 
leads to a power-law dependece : E l nea k oc E^ , with A ~ 0.57 
dAmatill2006tlAmati et alj|2002ll2009l) . 

There is also an additional phenomenological relation 
which has manifested itself in recent years, attracting at- 
tention: the possibility of fitting the observed spectra, in- 
tegrated over selected time intervals, by a blackbody plus 
power-law (BB+PL). There are cases in which the GRB spec- 
trum, again integrated over selected time intervals, is statisti- 
cally indistinguishable between fits with a Band function and 
a BB+PL; there are other cases in which the BB+PL model 
is even preferred over the Band function and finally there are 
sources f or which the data fulfill the Band functio n but not the 
BB+ PL (fRvda [200l l200l IRvde & Pe'eJ [20091: IRvde et alJ 
l20Toll201 lh . 

Generally, the additional blackbody componen t occurs in 
the ea rly part of the prompt emission (see e.g. IRvde et ail 
2010). There are bursts for which it is dominant with respect 
to the power-law, but opposite cases also exist, in which the 
power-law component is predominant with respect to the ther- 
mal one; th e relative intensity of the two components can vary 
with time dRvddl2004 [20051) . 

In some interesting cases the blackbody component shows 
a characteristic evolution in time, with both its observed tem- 



perature and flux having a dependence from the arrival time 
well des cribed by a broken power-law with index es aj, bj, af 
and b F (Ryde 2004, 20051 IRvde & Pe'er 200j). It is appro- 
priate to emphasize that the correspondence between differ- 
ent spectral models cannot, in general, be applied to the entire 
GRB light curve. There is then the necessity of identifying the 
specific time interval over which the different spectral features 
are identified, also possibly in order to re cover their physical 
origin (see e.g. Ruf fini et al.ll2010al.l201 lab . 

Specifically, for GRB 0803 19B and GRB 050904 both 
time-resolved and time-integrated BA T spe ctra, including the 
Proper-GRB (P-GRB, see below Sec. [3J) spectrum of GRB 
080319B, a re best fit by a simple power-law (see Sec. |4]and 
E\ see also iRacusinet al J 120081: IStamatikos et al.l l2009T and 
ICusumano et al. 1 120071) : alternative models, such as a Band 
function or a BB+PL, cannot be constrained by the BAT data. 
No concluding statements on the presence of a blackbody 
component in the P-GRB of GRB 0803 19B can be made due 
to the limited bandpass of the instruments (see Sec. [4). The 
analysis of Konus-WIND data from GRB 0803 19B has shown 
a best fi t with a Band functi on dRacusin et al.l2 008). For GRB 
050904 Sugita et al. (2009|) performed a joint spectral analy- 
sis among Swift BAT, Suzaku-WAM and Konus-WIND, find- 
ing a good fit with both a power-law with exponential cutoff 
and the Band function. A blackbody component has not been 
observed in either of these two sources. It is not clear at this 
stage if different conclusions could have been reached if these 
sources would have been observed by Fermi or AGILE. 

2.2. Theoretical models of GRB prompt emission 

Many different models have been developed to theoreti- 
cally explain the observational properties of GR Bs. One of 
the most quoted ones is the fireball model (see Piranl 120041 
for a review). An alternative one, originating in the elec- 
trodynamical processes a round a Kerr-Newmann black hole 
dDamour & Ruffini 19751) . is the fireshell model. 

The fireball model was first proposed b y | Cayallo & Reesl 
( fl978l) . iGoodmanl (fl986l) and iPaczvnskil d_1986j), who have 
shown that the release of a large quantity of 7-ray photons 
into a compact region can lead to an optically thick photon- 
lepton "fireball" through the production of e^ pairs. The term 
"fireball" refers to an opaque plasma whose initia l energy is 
significantly greater than its rest mass (lPiranlll999l) . 

The fireshell model also starts from an optically thick e^ 
plasma whose evolution has been followed in a sequence of 
states of thermal equilibrium, taking properly into account 
the ultrarelativistic expansion and the detailed com puta tion 
of the rate equation for t he e ± annihilation (see Sec. 13. II and 
iRuffini et al.lfi999ll2000l) . 

In the fireball model, the pr ompt emission, including 
the sh arp luminosity variations (Ramirez- Ruiz & Fenimorel 
2000), are due t o the prolonged and var i able activity of the 
"inner engine' ' dRees & Meszarosl[T994l IPiranl 12004). The 
conversion of the fireball energy to radiation occurs via 
shocks, either internal (whe n faster moving matter ov ertakes 
a slower moving shell, see Rees & Meszaroi 119941) or ex- 
ternal (when the moving matter is slowed down by the ex- 
ternal medium surrounding the burst, see iRees & Meszarosl 
119921) . Specifically, for GRB 0803 19B internal shocks 
have been considere d responsible for the prompt emission 
dRacusin et al.| [2008: see, however, Kumar & Naravan 2009) 
and external shocks are then c onsidered resp onsible for the 
afterglow (see also, howe ver. iKumari 119991 and references 
therein). For GRB 050904 ICusumano et al.l d2007l) proposed 
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that internal shocks are responsible f or both the prompt a nd 
the afterglow emission (see, however, Gendre et al. 2007 for 
a detailed description of limits and advantages of both the in- 
ternal and the external shock scenario). 

Much attention has been given to synchrotron emission 
from relativistic electrons, possibly accompanied by SSC 
emission to explain the observed GRB spectrum. These pro- 
cesses were found to be c onsistent with t he observational data 
of many GRBs (see e.g. lTavanilll996t iFrontera et al.ll2000t: 
Schaefer et al. 1998). However, several limitations have been 
evidenced in relation with t he low energy spec t ral slopes of 
time- i ntegrated spectra (see ICrider et al.l 119971 iPreece et alJ 
120021 IGhirlanda etaTI 120021 and IGhirlanda et al.l 120031; see 
also, ho wever, Daigne et al. 2009 ) and time-resolved spec- 
tra (see ICrider et al.l 119981 and TGhirlanda et all 12003b : ad- 
ditional limitations on SSC have also bee n poin ted out by 
IKumar & McMahonI d2008h and iPiran etaTI (120091) . 

In all the above considerations, the equations of motion 
of the fireball are evaluat ed under the ultrarelativist ic ap- 
proximation, leading to the Blandford & McKeel i 19761) self- 
similar power-law solution (see Sec. I3.lt . The maximum 
Lorentz factor of the fireball is estimated from the temporal 
occurrence of the peak of the optical emission, which is iden- 
tified with the peak of the forward external shock emission 
dMolinari et alj|2007t iRvkoff et al.ll2009l) in the thin shell ap- 
proximation "dSari & Piranlll999l) . It was also proposed to put 
an upper limit on the maximum fireball Lorentz gamma fac- 
tor from the upper limit on the intensity of a possible smooth 
background signal in the hard X-rays to soft gamma rays dur- 
ing the prompt emission, which is identified with the contri- 
bution of the forward ext ernal shock emission to the prompt 
phase dZou & Piranl2010T) . Another proposal was advanced to 
use compactness arguments within a scenario with two sep- 
arate emitting re gions for the MeV and the GeV emissions 
(IZou et alJl201 ll) . 

Partly alternative and/or complementary scenarios to the 
fireball model have been developed, e. g. the ones based 
on: quasi-thermal Comptonizatio n dGhisellini & Celoitil 
19991). Com p ton drag emission dZdziarski e t al.l 1 199 It 
Shemi 119941: lLazzati et al.l 120001) . synchrotron emissio n 
from a decayi ng magnetic fiel d dPe'er & Zhangl |2006), 
jitter radiation (Medvedev 2000), Compton scattering of 
synch r otron self-absorbed photo ns ([Panaitescu & Meszaios 
l2000t IStern & Poutanenl 120041), photospheric e mission 
dEichler & Levinsonl 120001: iMeszaros & Reesl 120001: 
Meszaios et al. 2002; Daigne & Mochkovitch 2002; Giannios 
20061: iRvde & Pe'er| j2009t lLazzati & Begelmanl 120101) . In 
particular, Ry de & Pe'erl (120091) pointed out that photo- 
spheric emission overcomes some of the difficulties of pure 
non-thermal emission models. 

In this paper we focus on the fireshell model. The charac- 
teristic parameters of the model are the total energy Ef ot , the 
baryon loading B and t he Circ umBurst Medium (CBM) den- 
sity n c b,„ (see Secs. l3.ll and l3.4t . The Lorentz 7 factor, directly 
linked to B and E* ot , is explicitly computed from the descrip- 
tion of all the phases starting from the moment of gravitational 
collapse. The fireshell radial coordinate is explicitly evaluated 
as a function of the laboratory time, the comoving time and 
the arrival time at the detector (see Sec. I3.lt . The optically 
thick e ± plasma, endowed with baryon loading, is followed 
all the way to transparency. The collision with the CBM of 
the emerging relativistic expanding shell of baryons gives rise 
to the extended afterglow, which comprises both the prompt 



emis sion and the traditional decaying afterglow phases (see 
Sec. I3.lt . Relativistic effects are taken into account in the 
computation of the equations of mo tion of the shell and of 
the EO uiTemporal Surfaces (EQTS, iBianco & Ruffinill2004l 
2005b, see Sec. 13 . 3b . Taking into proper account these rel- 
ativistic effects, it is possible to deduce the spectrum of the 
collision process between the baryons and the CBM in the 
comoving frame of the shell. For simplicity it was initially 
assumed that such collisions occur in a fully radiative regime 
and give rise to a p ure thermal spectrum in the comoving 
frame (see Sec. 13. 3t . The observed spectrum is then ob- 
tained as a double convolution of thousands of thermal spec- 
tra, each one with a different temperature and weighted by 
the appropriate Lorentz and Doppler factors, following the 
solution of the equations of motion of the fireshell, both over 
the EQTS and over th e observation time dRuffini et al.ll2004l 
IBernardini et"al1l200M 

As we recalled in the Introduction, this ansatz of a pure 
black body spectrum in the fireshell comoving frame, in 
spite of its simplicity, led to a successful interpretation of 
many different sourc es with an E ob e s ak ins ide the instrumental 
bandpass (see e.g. Bernardini et al. 2005; Ruffini et aLl [2006 
Dainotti et al. 2007; Bernardini et al. 2007; Caito et al. 2009 



2010; de Barr os et all 1201 ll) . We were then quite confident 



that we could obtain any observed power-law by an adequate 
convolution of thermal spectra duly weighted by the relativis- 
tic Doppler factors. The impossibility of obtaining the correct 
power-law indexes observed in both the sources treated in this 
paper, which have an E° b e s ak above the instrumental bandpass, 
using solely thermal spectra in the fireshell comoving frame, 
was quite unexpected. This impossibility stems from the fact 
that neither the distribution of the temperature in the comov- 
ing frame, nor the Doppler factors used in the thousands of 
convolution processes, can be arbitrarily given. Both of them 
are constrained by the equations of motion of the fireshell 
and by the consequent release of kinetic energy in the colli- 
sion. A priori, a convolution of thermal spectra with arbitrary 
temperatures and Doppler factors can always fit any observed 
power-law. There is, however, no way to fit by convolutions 
of thermal spectra the observed power-laws consistently with 
the fireshell equations of motion in these two sources with 
E° b e s ak above the instrumental bandpass. 

What is particularly interesting, however, is the fact that this 
difficulty can be overcome simply by adding an extra power- 
law component to the pure black body spectrum. As we will 
see, the application of the fireshell model to GRB 0803 19B 
and GRB 050904 leads in fact to the introduction of an ad- 
ditional phenomenological parameter a, which characterizes 
the departure of the slope of the low energy part ofthe co- 
moving spectrum from a pure thermal one (see Sec. 13.3.21 . 

This new result is consistent also with all the previously an- 
alyzed GRBs, and reasons for this are given in Fig. Q] The 
success of this approach is not trivial, since there is no direct 
analytic relation between the index of the power-law intro- 
duced in the spectrum in the fireshell comoving frame and the 
observed one. It is quite significant that the introduction of 
a single power-law makes the fireshell model consistent with 
all observed GRB spectra. 

The physical explanation for a has not yet been found. 
Analogously, no physical explanation has yet been found for 

6 Typically, a resolution of ~ 5 X 10 thermal spectra for each second of 
observation has been used. 



Analysis of GRB 0803 19B and GRB 050904 within the fireshell model: evidence for a broader SED 



5 



the above described (Sec. 12.1b phenomenological parame- 
ters of the Band function and of the Ama t i relat i on, an d in 
this sense for the ones described by iRvdd (120041 120051) and 
Ryde & Pe'er (2009) as well. All these phenomenological 
parameters lead to a better quantitative description of GRBs; 
they are the object of active theoretical studies and are an im- 
portant step toward reaching the future identification of the 
underlying physical processes characterizing the GRBs. 

We now proceed to a detailed description of the fireshell 
model (Sec. |3]l. 

3. THE FIRESHELL MODEL: THE GRB LUMINOSITY AND 
SPECTRUM 

The fireshell model, avoiding a piecewise fit of the observed 
GRB data, proposes a unified picture starting from the ini- 
tial process of gravitational collapse to a black hole through 
all successive stages. It gives a theoretical treatment of each 
stage, identifying the parameters intrinsic to the source, essen- 
tial to describing the evolution of the system, as well as its in- 
teraction with the CBM. The corresponding equations of mo- 
tion are treated accordingly. Regimes of Lorentz gamma fac- 
tors in the range 100-1000 are encountered, implying the ne- 
cessity of a fully relativistic approach. The model is intrinsi- 
cally endowed with highly nonlinear effects: each prediction 
of the theoretical model at a given instant of the observation 
time is influenced by the entire history of the source. Conse- 
quently, any solution of the model in agreement with the ob- 
servations necessarily implies a high level of self-consistency. 

3.1. The canonical GRB 

Within the fireshell model, all GRBs originate from an op- 
tically thick electron-positron plasma in thermal equilibrium, 

having total energy Ef n , and formed in the gravitational col- 
lapse to a black hole (Ruffi niet alJ l2010bh . The condition 
of thermal equ i libriu m assumed in our model and proved by 
lAksenov et al.l (120071) distingui shes our model from alterna- 
tive approaches (e.g. the one bv lCavallo & ReesTT9 78). where 
the total annihilation of the e ± plasma was assumed, leading 
to a vast release of energy pushing on the CBM (the concept 
of a "fireball"). In our case the annihilation of the e pairs 
occurs gradually and is confined within an expanding shell: 
the "fireshell". The rate equation for the e ± pairs and their 
dynamics (the pair-electr omagnetic pu l se or PEM pulse for 
short) has been given by Ru ffini et alJ (119991) . This plasma 
engulfs the baryonic material left over in the process of 
gravitational collapse having mass Mb, still keeping thermal 
equilibrium between electrons, positrons and baryons. The 
baryon loading is measured by the dimensionless parameter 

B = M b c 2 /e£, . It was shown (see IRuffini et alJl2000l) that no 
relativistic expansion of the plasma can be found for B > 10~ 2 . 
The fireshell is still optically thick and self-accelerates to 
ultrarelativistic velocities (the pa ir-electromagnetic -b aryonic 
pulse or PEMB pulse for short, IRuffini et all 120001) . Then 
the fireshell become s transparent and the P-GRB is emitted 
(IRuffini et al.ll2001al) . The amount of energy radiated in the 
P-GRB is only a fraction of E* ot . The remaining energy is 
stored in the kinetic energy of the optically thin baryonic and 
leptonic matter fireshell. The final Lorentz 7 factor at trans- 
parency, 70, can vary in a vast range between 10 2 and 10 3 as 

a function of Ef% and B (IRuffini et alj|2000h . 

After transparency, the remaining accelerated baryonic 
matter still expands ballistically and starts to slow down by 
the collisions with the CBM, having average density n c i, m . 



During this phase, t he extended afterglow emission occurs 
( IRuffini et alj|2001al) . In common with the majority of exist- 
ing models, we describe the motion of the baryons as an ex- 
panding thin shell enforcing energy and momentum conserva- 
tion in the collision with th e CBM. The condit ion of a fully ra- 
diative regime is assumed (Ruffini et al. 2003). It is appropri- 
ate to recall a further difference between our treatment and the 
ones in the current literature. The complete analytic solution 
of the eq uations of motion of the baryon ic shell has been de- 
veloped ( Bianco & Ruffini 2004, 2005b), while i n the current 
literature usually the Blandford & McKed (Il976|) self-similar 
solution has been uncriti cally adopted (e . g. | Meszaros et al.| 

' | 1998t IWaxmad Il997t iRees & Meszaros 

1999tlPanaitescu & Mesza roslll998U Piran 
1999; Gruzinov & Waxman 1999; van Paradiis et al.l l2000t 
Meszaros 2002). The similarities and differences between 
the two approaches hav e been explicitly pointed out in 
iBianco & Ruffinil (l2005al) . 

From this general approach, a canonical GRB bolomet- 
ric light curve is defined which is composed of two differ- 
ent parts: the P-GRB and the extended afterglow. The rela- 
tive energetics of these two components, as well as the ob- 
served temporal separation between the corresponding peaks, 
is a function of the above three parameters E e tot , B, and « c /, m ; 
the first two parameters are inherent to the accelerator char- 
acterizing the GRB, i.e., the optically thick phase, while the 
third one is inherent to the GRB surrounding environment 
which gives rise to the extended afterglow. What is usually 
called the GRB "prompt emission" in the literature is ac- 
tually composed of both the P-GRB and the initial part of 
the exte nded afterglow encom passing its peak. As we pro- 
posed in Ruffini et al.l (12001a). both the so-called "short" and 
"long" GRBs fit into this canonical GRB scenario. In par- 
ticular, for baryon loading B < 10~ 5 , the P-GRB component 
is energetically dominant over the extended afterglow. In the 
limit B — >• it gives rise to a "genuine short" GRB. Other- 
wise, when 3.0 x 10~ 4 < B < 10~ 2 , the kinetic energy of the 
baryonic and leptonic matter, and consequently the extended 
afterglow emission, predominates with respect to the P-GRB, 
giving rise to the "long" GRBs or the "disguised short" GRBs 
depending on the avera ge CBM densi ty and the astrophysical 
scenario (IRuffini et alJpOOlal 120021; iBernardini et al.l 12007b 
ICaito et al.l2009Ll2010tlde Barros et al.l201 ll) . Since the "crit- 
ical" value of B (i.e. the value of B for which both the P-GRB 
and the extended afterglow have the same energy) is a slowly 

varying function of E e m , for 10" 5 < B < 3.0 x 10" 4 the ratio of 
the total energies of the P-GRB and of the extended afterglow 
is also a function of Ef ol (Ruffini et al. 2009). 

If one goes to the observational properties of this model of 
a relativistic expanding shell, a crucial concept has been the 
introduction of the EQTS. In this topic, also, our model is 
distinguished from those in the literature for deriving analytic 
expressions for the EQTS from the analytic solutions of the 
equations of motion (Bianco & Ruffini 2005b). 

The observed temporal variability of the extended afterglow 
is produced in our m odel by the interac tion of the fire shel l 
with CBM "clumps" (IRuffini et alJl2002l see also Sec. EH- 
The issue of time variability in GRB li ght curves has bee n 
longly debated. Sever al authors, e.g. Zhang e t al.l (12006b : 
iNakar & Granotl (120071) . found that CBM inhomogeneities are 
not able to produce the short-timescale variability in GRB 
prompt and afterglow emission light curves, mainly because 
photons emitted at the same instant of time from different 
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parts of the emitting regions have different arrival times: this 
tends to smoothen the light curve significantly. 

Within the fireshell model, also on this point there are some 
differences with the other models. It is emphasized that, from 
the correct computations of the equations of motion of the 
shell and of the Lorentz 7 factor, the short time scale vari- 
ability of GRB light curves occurs in regimes with the largest 
values of the Lorentz gamma factor, when the total visible 
area of the emission region is very small and the "dispersion" 
in arrival time of the luminosity peaks is negligible. There- 
fore, under this condition the short- timesc ale variability of 
GRB light curves can be produced by inhomoge neiti es in the 
CBM, as fou nd also bv lDermer & Mitmanl ( 1 19991) and lDermerl 
(120061 120081) . The application of the fireshell model leads to 
a direct evaluation of the filamentary and clumpy structure of 
the CBM, which was already predicted in pioneering works 
by Enrico Fermi in the t heoretical study of Interstellar Mat- 
ter (ISM) in our galaxy dFermill 19491 119541) . and is much on 
line with the knowledge obtained f rom various studies of the 
ISM in galaxies (see, for exa mple, iKim et alJll998t iLockmanl 
20021 iFinkelstein et al.ll2009l) . 

Clearly, the fact that the short time scale variability is ob- 
served only in the prompt emission, while the X-ray after- 
glow light curves are usually smooth, does not mean that 
the CBM is inhomogeneous only up to a given radius, be- 
yond which it becomes homogeneous. Inhomogeneities are 
everywhere in the CBM, but beyond a given distance from 
the source, corresponding approximately to the end of the 
prompt emission phase, they do not produce observable ef- 
fects on the light curve since they are indeed smeared out 
by the curvature effect and the relativistic effects between 
the time in the fireshell comoving frame and the photon ar- 
rival time at the detector (Ruffini et al. 2002, 2006). The early 
X-ray afterglow originates in fact from the same kind of in- 
teraction of the fireshell with a clumpy CBM. The absence 
of spiky emission is simply a consequence of these effects. 
In other words, it is true that the fireshell model predicts 
that the same clumps at lar ger radii would produce longer 
spikes dRuffini et alj 1200 lbh . However, at such large radii 
where the effect would be measurable, the s mearing dom- 
inates and p revents the effect to be observed (Ruffi ni et al.l 
2OO2II2OOI). Vice versa, the prompt emission, where the ef- 
fect of the CBM clumps are observable, occurs encompass- 
ing too limited a range of radial distances to make this ef- 
fect noticeable in all sources. The only exception may oc- 
cur in the case of isolated high density clumps along the line 
of sight at late time of emission where indeed this effect is 
observable (iBianco et al.ll2Q06t iBernardini et al] r2008. 2009; 
iRuffini et al.ll2009Hlzzo et al.ll2Q10t) . We are currently verify- 
ing if this aspect of the fireshell model may also explain some 
specific properties of the X- ray flares recently ev idenced in 
the afterglow phase (see e.g. Margutti e t al.ll201 lallbl and ref- 
erences therein). The drop in energy of 2-4 orders of mag- 
nitude from the prompt 7-ray phase to the early X-ray af- 
terglow is perfectly in line with the decrease of the Lorentz 
gamma factor in the expansion of the fireshell. The fact that 
the prompt emission in 7-rays stops at some time is related 
to the well known hard-to-soft evolution of the emission pro- 
cess, wh ich i s perfectly explained within the fireshell model 
(see Sec. 13.31) and it is related to the solution of the equations 
of motion of the system, implying the decrease of the Lorentz 
gamma factor as well as the amount of energy release in the 
collision between the fireshell baryonic m atter and the CBM 
(IRuffini et al.ll2004t fBernardi ni et al.ll2005l) . 



3.2. The P-GRB spectral properties 

The spec tr um a t transparency has been given in 
IRuffini et al.l (|2000) with a temperature computed con- 
sistently with the local thermodynamics of the e ± -baryon 
plasma and the Lor entz gamma factor. Details are given in 
IRuffini et al.l (|2009). It is appropriate to stress that in the 
emission of the P-GRB there are two different contributions: 
one corresponding to the emission of the photons due to the 
reach of the transparency, and the second originating from 
the interaction of the protons and electrons with the CBM. A 
spectral energy distribution with a thermal component and a 
non-thermal one should be expected to occur. 

3.3. The extended afterglow spectral properties 

The majority of work in the current literature has ad- 
dressed the analysis of the prompt emission as originating 
from various co mbination s of synchrotron and inverse Comp- 
ton processes dPiranl 120041) . It appears clear, however, that 
this interpretation is not satisfactory (see Sec, [p see also 
Ghirla nda et al.ll2003b iKumar & McMahodl2W8HPiran et alJ 
2009). Furthermore, in the description of an ultrarelativistic 
collision between protons and electrons and the CBM new 
collective processes of ultrarelativistic plasma physics occur, 
not yet fully explored and un derstood (e.g. Weibel instability, 
see IMedvedev & Loeb|[T999l) . Most promisin g resu l ts alon g 
this line have been already ob tained by ISpitkovskyl (|2008a); 
IMedvedev & SpitkovskyI (120091) . 

Without waiting for the developments of these investiga- 
tions, we have adopted a very pragmatic approach in the 
fireshell model by making full use of the knowledge of the 
equations of motion, of all the EQTS formulations as well 
as of the correct relativistic transformations between the co- 
moving frame of the fireshell and the observer frame. In this 
respect, we have adopted a fundamental procedure: to make 
an ansatz on the spectral properties of emission of the col- 
lisions between the baryons and the CBM in the comoving 
frame, and then evaluate all the observational properties in 
the observer frame. In order to take into proper account the 
filamentary, clumpy and porosity structure of the CBM, we 
have introduced an additional parameter 1Z, which describes 
the fireshell surface filling factor. It is defined as the ratio be- 
tween the effective e mitting area of the fireshell A e ff and its 
total visible area A vis (IRuffini et al.ll2002l 120041 120051) . 

It must be emphasized that the fact that only a fraction 1Z of 
the shell surface is emitting does not mean that only a fraction 
1Z of the total shell energy is emitted. We must in fact dis- 
tinguish between an instantaneous interaction of the fireshell 
with a single filament and its overall interaction with all the 
filaments of the entire cloud giving rise to the spiky struc- 
ture of the light curve. This global interaction is clearly the 
superposition of randomly distributed instantaneous events. 
The different filaments inside the cloud interact with different 
parts of the fireshell and the entire cloud reduces the kinetic 
energy of the entire fireshell. The key point is that, during the 
prompt emission, the cloud, typically with a mass of the order 
of 10~ 8 - 10 -11 solar masses, covers the entire visible area of 
the fireshell (typically with a radius between 10 13 cm and 10 15 
cm for the sources analyzed in the present paper, see Tables Q~| 
and|2). Consequently, at any given instant of time, each fila- 
ment of the cloud covers only a small fraction of the fireshell 
surface. However, when we integrate over the cloud crossing 
time, the coverage of the cloud as a whole is equal to unity. 

As a first ansatz, we have assumed that the extended af- 
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terglow radiation has a thermal spectru m in the comoving 
frame of the fireshell dRuffini et alj2004l) . The observed GRB 
spectrum is given by the convolution of hundreds of thermal 
spectra with different temperatures and different Lorentz and 
Doppler fact ors. Such a convolution is to be performed over 
the EQTSs dBianco & Ruffini 2004, 2005b), which are the 
surfaces of constant arrival tim e of the photons at the d etector, 
and over the observation time (Bemardini et al. 2005). 

The hard-to-soft transition of GRB time-integrated and 
time-re solved spectra, th at was first observed in BATSE 
GRBs (ICrider et al.l fl997). within the fireshell model comes 
out naturally from: 1) the evolution of the comoving tem- 
perature; 2) th e decrease of the bulk 7 factor; 3) the cur- 
vature effect dBianco & Ruffinil l2004t iRuffini et alj 12004 
iBernardini et al.ll2005l) . 

Within the fireshell model, the extended afterglow luminos- 
ity at the detector arrival time t% per unit of solid angle dCt and 
in the energy band [v\ , u 2 ] is given by (Ruffini et al. 2004) 



EQTS 



^vcos0 A 4 ^W(iy u is 2 ,T arr )dX, (1) 
47r at" 



where Ae = AEj„,/V is the emitted energy density released 
in the interaction of the accelerated baryons with the CBM 

measured in the comoving frame, A = {^[l - {v / c)cosd]} 1 
is the Doppler factor, ^(y\,v 2 ,T arr ) is an "effective weight" 
required to evaluate only the contributions in the energy band 
[v\,vi\, dY, is the surface element of the EQTS at detector 
arrival time t d a on which the integration is performed and T arr 
is the observed temperature of the radiation emitted from t/E. 
The "effective weight" ^N{v\,v 2 ,T arr ) is defined as the ratio 
between the energy density emitted in a given energy band 
[vi , v 2 \ and the bolometric energy density: 



W(vi,v 2 ,T arr )-- 



ii 



00 / clN^ 
\dVde 



ede 



(2) 



where is the number density of the photons per unit of 
energy in the comoving frame of the fireshell. 

3.3.1. Thermal case 

With the assumption of a thermal spectrum in the comoving 
frame of the fireshell 



aW 7 _ / 8tt 
dVde ~ \ 



(3) 



(h is the Planck constant, c is the speed of light and kg is the 
Boltzmann constant) we have 



W(y u v 2 ,T arr )-- 



aT 4 



(4) 



where a is the radiation constant and T is the temperature in 
the comoving frame. 

In general, the temperature in the comoving frame can be 
evaluated starting from the following relation: 



AE„ 



= Trr 2 cll I ( -Jh- I tde, 



(5) 



At "' ~"J \dVde, 
where At is the time interval in which the energy AEi nt is 



E 

o 

a3 
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Figure 2. Theoretically simulated instantaneous spectra obtained by assum- 
ing £™_ = 1.0 X 10 54 erg, B = 2.5 X 1(T 3 , n cbm = 1 part cm" 3 and different 
values of the index a. The curve with a = 0.0 corresponds to the pure thermal 
spectrum case. 

developed. By inserting Eq. © into Eq. (0 we obtain 

T= ^vrr^Ar) ' (6) 
with a the Stefan-Boltzmann constant. 

3.3.2. The ' 'modified" thermal spectrum 

The new SED for the radiation emitted in the comoving 
frame of the fireshell is a "modified" thermal spectrum: a 
spectrum characterized by a different asymptotic power-law 
index in the low energy region with respect to the thermal 
one. This index is represented by a free parameter a, so that 
the pure thermal spectrum corresponds to the case a = 0: 



dN-y 
dVde 



87T 



e 



exp 



(7) 



a is a phenomenological parameter defined in the comoving 
frame of the fireshell. This phenomenological approach can 
be relevant in identifying the true physical mechanisms oc- 
curring in the collisions in the comoving frame, and uniquely 
separates them from the relativistic contributions coming 
from relativistic transformations and convolutions over the 
EQTS leading to the observed spectrum. 

By using the Eq. and introducing the variable y = 
e/(kgT), we obtain the following expression for the "effective 
weight": 



W{v\,v 2 ,T arr )-- 



r(4 + a)L/ 4+Q (l) 



(8) 



where T(z) = / °°f : '<? 'dt is the Gamma function and 
Li n (z) = z k /k n is Jonquiere's function. 

In analogy with the thermal case, we can define an "effec- 
tive temperature" for the "modified" thermal spectrum; by in- 
serting Eq. (O in Eq. (0 we obtain: 



T = 



f AE it 



3, .2 



he 



1/4 



\ At J (4irr 2 )2irllk 4 B r(4 + a)Li4 +a (l)_ 
It can be easily seen that, for a = 0, we obtain Eq. (O. 



(9) 
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In Fig. |2] are shown several theoretically predicted instanta- 
neous spectra characterized by the same temperature and dif- 
ferent values of the free parameter a. It can be seen that the 
main effect of varying the value of a is a change in the low en- 
ergy slope of the spectral energy distribution. In particular, by 
decreasing a the low energy emission increases. Around the 
peak energy the spectrum is instead only weakly dependent 
on the value of a. 

In the present paper we test the comoving modified thermal 
spectrum by comparing the numerical simulations with the 
observed prompt emission spectra and light curves of highly 
energetic GRBs. In particular we present, as specific exam- 
ples, the analysis of the observational BAT data (in the 15- 
150 keV bandpass) of G RB 0803 19B (E iso = 1.34 x 10 54 erg, 
iGolenetskii etai]|2008l) and GRB 050904 (E iso = 1.04^ x 
10 54 erg. Sugit a"eTaT]l2009l) in Secs.|4]and|5] respectively. 

3.4. The numerical simulation of GRB light curves and 
spectra 

To best reproduce the observational data within the fireshell 
model, we need to determine the following five parameters: 

Ef ot , B, a, n cbm and TZ. 

The procedure assumes a specific value of E e tot and B. It 

is clear that E e tot has to be larger or equal to the observed 

isotropic equivalent energy Ej SO of the GRB. Ef ot can be ac- 
tually quite larger than Ej SO since, in many sources, we are 
limited by the threshold and the bandpass of the detectors. 
The value of B is determined by the ratio between the ener- 
getics of the P-GRB and of the extended afterglow, as well 
as by the time separation between the correspo nding peaks 
(Ruffi ni et al.ll2001all2008l:lAksenov et alJl20Toh . 

The determination of the three remaining parameters de- 
pends on the detailed "fitting" of the shape of the extended af- 
terglow light curves and spectra. In particular, the parameter 
TZ determines the effective temperature in the comoving frame 
and the corresponding peak energy of the spectrum, a deter- 
mines the low energy slope of the comoving spectrum and 
n c bm determines the temporal behavior of the light curve. It is 
found that the CBM is typically formed of "clumps" of width 
~ 10 15 - 10 16 cm and density contrast 10" 1 <Sn/n< 10. Par- 
ticularly important is the determination of the average value 
of n c b m . Values of the order of 0. 1-10 particles/cm 3 have been 
found for GRBs exploding inside star forming region galax- 
ies, while values of the order of 10~ 3 particles/cm 3 have been 
found for GRBs explodi ng in galactic ha l os (i.e. the "dis- 
guise d" GRBs. see iBernardini et al.l l2007t ICaito et al.ll2009l 
2010: lde Barros et al.H201 ilT ~ 

Of course, "fitting" a GRB within the fireshell model is 
much more complex than simply fitting the N(E) spectrum 
with phenomenological analytic formulas for a finite temporal 
range of the data. It is a consistent picture, which has to "fit" 
the intrinsic parameters of the source, as well as its spectrum 
and its light curve temporal structure. Concerning the theoret- 
ical spectrum to be compared with the observational data, it 
is obtained by an averaging procedure of instantaneous spec- 
tra. In turn, each instantaneous spectrum is linked to the fit of 
the observed multiband light curves in the chosen time inter- 
val. Therefore, both the "fit" of the spectrum and of the ob- 
served multiband light curves have to be performed together 
and jointly optimized. Moreover, the parameters used in the 
numerical simulations are not independent. In fact, they have 
to be computed self-consistently through the entire dynami- 



cal evolution of the system and not separately at each time 
step. For each spike in the light curve the parameters of the 
corresponding CBM clumps must be computed, taking into 
proper account all the thousands of convolutions of comoving 
spectra over each EQTS leading to the observed spectrum. It 
is clear then that since the EQTS encompass emission pro- 
cesses occurring at different comoving times, weighted by 
their Lorentz and Doppler factors, the "fitting" of a single 
spike of the light curve is not only a function of the prop- 
erties of the specific CBM clump but of the entire previous 
history of the source. Any step of the "fitting" process affects 
the entire following evolution and, viceversa, at any step a 
"fit" must be made consistently with all the previous and sub- 
sequent history: due to the non linearity of the system and to 
the EQTS, any change in the "fit" produces observable effects 
up to a much later time. This implies that the "fitting" process 
cannot proceed for successive temporal steps: the complete 
analysis must be applied to the entire GRB as a whole, to 
avoid possible systematic error propagation from a temporal 
step to the following ones. This leads to an extremely com- 
plex trial and error procedure in the fitting of the data in which 
the uniqueness of the parameters defining the source are fur- 
ther and further narrowed down. Of course, we cannot expect 
the latest parts of the fit to be very accurate, since some of the 
basic hypotheses on the eq uations of motion, a nd the possible 
fragmentation of the shell dDainotti e t al. 2007), can affect the 
fitting procedure. 

4. GRB 0803 19B 

We analyzed the GRB 0803 19B prompt emission light 
curve and spectrum observed by BAT within the fireshell 
model. As we already mentioned in the introduction, for GRB 
0803 19B we have E° b e s ak = 675 ± 22 keV, although there is a 
hard-to-soft spectral evolution going from with E"^ e s ak = (748 ± 
26) ke V at 22 s after the BAT trigger to Ef e s ak = (528 ± 28) ke V 
at 24 s after the BAT trigger, and the BAT spectral index re- 
duces from ~ 1 .0 to ~ 2.1 at 53 seconds after the BAT trigger 
dRacusin et al.ll2008l) . 

Several authors found some evidence of the possibility to 
separate the prompt emission of this source into two main 
episodes, parti t ioned at about 28 s after the BAT trigger time. 
Marg utti et al.l (120081) analyzed the variability time-scale t var 
of the 7-ray prompt emission, finding that the first part of the 
light curve (up to ~ 28 s) is dominated by t var ~ 0. 1 s, while 
the last part s hows a much longer ch aracteristic time-scale 
(t var ~ 0.7 s). IStamatikos et al.l (120091) found that the arrival 
offset between the Swift-BAT 15-25 keV and 50-100 keV en- 
ergy band (7-ray spectral lag) is maximum at t > 28 s and it 
appears to be anti-correlated with the arrival offset between 
prompt 15-350 keV 7-rays and the optical emission observed 
by TORTORA (optical/7-ray spectral lag), maximum at t < 
28 s. 

Concerning the first episode, we identify the first 7 s of 
emission (from -5 s up to 2 s after the BAT trigger time) 
with the P-GRB; the theoretically estimated total isotropic en- 
ergy emitted in the P-GRB and the observed temperature are 
E^ CRB = 1.85 x 10 52 erg and T°^ RB - 16 keV respectively. 
There are three main reasons that supports this interpretation: 

1. First, we performed the analysis of BAT spectra in- 
tegrated over sub-intervals of time encompassing the 
whole prompt emission by using the standard FTOOLS 
package (Heasoft, version 6.10). We found that all 
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these spectra are best modeled with power-laws and a 
discontinuity in the hard-to-soft evolution came out a 
few seconds aft er the BAT trigger tim e, as shown in 
Fig. [3] (see also Stamat ikos et al] l2009S). In particular, 
there is a clear soft-to-hard evolution up to ~ 1 s after 
the BAT trigger time, while a typical hard-to-soft tran- 
sition starts at about 8 s after BAT trigger time. It is 
difficult to evaluate at what time exactly the discontinu- 
ity occurs: in fact, in the region between ~ 1 s and ~ 8 
s, the photon index first appears to reach an asymptotic 
value of ~ 0.76, then further decreases up to ~ 0.7: this 
behaviour could be due to the partial superimposition of 
the contributions of both the P-GRB and the extended 
afterglow. 

2. The second reason suggesting the interpretation of the 
first seconds of emission as the P-GRB is that the opti- 
cal emission starts at t ~ 9 s after the BAT trigger time 
(see Fig.|4]l: in fact, within the fireshell model the opti- 
cal radiation is expected in the extended afterglow, but 
not in the P-GRB. 

3. It is important to point out that, besides the two above 
observational considerations, the identification of the 
first 7 s of emission as the P-GRB is the only interpre- 
tation that allow us to constrain the values of Ef ot and 
B in such a way as to make a consistent fit of both the 
P-GRB and the extended afterglow. 

Of course, none of these three arguments separately would 
give support to our interpretation: e.g., the optical component 
could have also been present before ~ 9 s but unobserved due 
to the instrumental constraints. However, the redundant oc- 
currence of all three arguments implies that our interpretation 
of the first 7 s as the P-GRB component is fully compatible 
both with our theoretical framework and with the observed 
data. The P-GRB spectrum is expected to be co mpos ed of 
a thermal plus a non-thermal component (see Sec. 13.21 1. We 
found that the BAT spectrum of the first 7 s of the prompt 
emission is best fit by a power-law with photon index 7=0.84 
± 0.04, with a chi square value of \ 2 = 48.26 for 60 degrees 
of freedom. The presence of the expected additional thermal 
component cannot be constrained by the data (C. Guidorzi, 
private communication). A possibility is that thermal flux is 
much lower than the non-thermal one and then it is negligi- 
ble. Another alternative and/or complementary possibility is 
that the thermal component has been missed due to the limited 
bandpass of the instruments. It can be only matter of specula- 
tion if observations by Fermi and AGILE would have led to a 
different conclusion. 

The remaining part of the first episode (from ~ 7 s up to 
about 28 s) is interpreted as the peak of the extended af- 
terglow, whose temporal variability is produced by the in- 
teraction with the CBM. The numerical simulation that best 
reproduces the light curve (Fig. |5) and the time-integrated 
spectrum (Fig. |7]i of this first episode (3s <t^< 28s) is ob- 
tained with the following parameters: Ef ot = 1.32 x 10 54 erg, 
B = 2.3 x 10~ 3 and a = -1.8; the Lorentz gamma factor at 
the transparency point, occurring at ro = 2.8 x 10 14 cm, is 
7o = 428. We consider an average number density (n e / jm ) ~ 6 
particles cm" 3 and K=3.5x 10" 10 . The structure of the CBM 
adopted is presented in Fig. [6] and the adopted density con- 
trast with respect to the average density is reported in Tab. Q] 
The distribution of the CBM is just an approximation of the 
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Figure 3. Time evolution of the photon index for power-law fit to 15-150 keV 
Swift BAT spectra integrated over time intervals of 3 s. A clear soft-to-hard 
evolution up to ~ 1 s after the BAT trigger time, while a typical hard-to-soft 
transition starts at about 8 s after the BAT trigger time. In the region between 
~ 1 s and ~ 8 s, the photon index first appears to reach an asymptotic value 
of ~ 0.76, then further decreases up to ~ 0.7. 
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Figure 4. Comparison of the prompt optical light curve (TORTORA, blue 
points; RAPTOR-Q, green points; Pi of the Sky, blue points) and the hard 
X-ray to 7-ray light curve (IBAS, grey points; BAT, black points) of GRB 
0803 19B; the pink dotted line marks the begin of the optical emission de- 
tected by TORTORA (t ~ 9 s after the BAT trigger time), whose onset is 
very rapid (~ t 3 08 ). Reproduced from Wozniak et al. (2009) with kind per- 
mission of PR. Wozniak and of AAS. 

real one, where the CBM density shows some smooth fluctu- 
ations around its trend during the fireshell evolution. Never- 
theless, it is sufficient to account for the observed variability 
in the luminosity. We must note that there is a sharp and short 
spi ke in the light curve 12. 4 s after the BAT trigger time (see 
e.g. lStamatikos e t al. 2009) which we are unable to reproduce 
within our model, based on a spherically symmetric approxi- 
mate dynamics. For this spike to be interpreted, a fully three- 
dimensional description of the CBM is needed. However, we 
expect that this more detailed description will not modify the 
overall dynamics of the system. In fact, the fluence of this 
spike is ~ 4.9% of the fluence observed in the first episode 
(between 2 s and 28 s) and ~ 2.5% of the fluence of the entire 
prompt emission (between 2 s and 57 s). Therefore, the error 
introduced by the omission of this spike from the numerical 
simulation is much smaller than the difference between the 
co-moving pure thermal spectrum and the modified one dis- 
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Figure 5. Theoretically simulated light curve of GRB 0803 1 9B prompt emis- 
sion in the 15-150 keV energy band (black solid curve) is compared with the 
data observed by BAT (red points); the P-GRB is marked with a magenta cir- 
cle. The vertical dotted line marks the begin of the second part of the prompt 
emission (t ~ 28 s). The labels "a", "b", "c", "d", "e", "f", "g" and "h" 
identify the peaks (see Fig.|6]and Tab.[T}. 
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Figure 6. Structure of the CBM adopted for GRB 0803 19B. The labels "a", 
"b", "c", "d", "e", "f", "g" and "h" indicate the values corresponding to the 
peaks in the BAT light curve (see Fig.[5]and Tab.[T}. 
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Table 1 

Properties of the CBM structure adopted for GRB 0803 19B: distance from 
the center of the explosion, thickness, normalized density and mass of the 
clumps; for each distance the transverse dimension of the visible area is also 
reported. 
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Figure 7. Theoretically simulated spectra of GRB 0803 19B integrated over 
the time interval 3s < t% < 28s with a = — 1.8 (black solid line) and a = 0.0 
(pure thermal case, green dashed line) are compared with the data observed 
by BAT (red points). It can be seen that with the "modified" thermal spectrum 
we can correctly reproduce the observed spectrum, contrary to what happens 
with the pure thermal spectrum. In the residual plot the pure thermal case is 
omitted. 



cussed in this paper. We can then omit the spike from our 
analysis, without hampering qualitatively the conclusions of 
our paper. 

With the above described set of parameters it is possible to 
interpret also successfully the spectra integrated over smaller 
intervals of time. Fig[8]shows, as an example, the spectrum for 
3s < fa < 13s: it can be seen that with the modified thermal 
spectrum we can correctly reproduce also this spectrum; on 
the contrary, by assuming a comoving thermal spectrum there 
are several discrepancies between the theoretical prediction 
and the observational data, especially at the lower energies. 
This is an important check to be made each time. In fact, 
changing the spectrum integration time means changing the 
number of different co-moving spectra which are convolved 



to get the observed one. The fact that the model is able to 
reproduce the observed spectrum regardless of the time scale 
over which it is integrated is therefore a clear support of the 
correctness of the assumed co-moving spectral shape. 

To have an estimate of the sensitivity of the determina- 
tion of these parameters of the model, we can proceed as 

follo ws. We fix Ef„, to the observed value of 1.32 x 10 54 
erg dGolenetskii et al.l 120081) . From the observational data 
we have that the fluence of the first 7 seconds, which cor- 
respond to the fluence of the P-GRB, f P - G RB, is 2.19 x 10" 6 < 
fp-CRB ^$ 2.29 x 10~ 6 erg/cm 2 . This fixes a range of values for 
B: 2.19 x 10" 3 < B < 2.33 x 10" 3 . Correspondingly, we must 
have 1.5 x 10" 10 < K < 6.0 x 10" 10 and 4.1 < (n chm ) < 8.2 
particles/cm 3 to reproduce the observed light curves and spec- 
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Figure 8. Theoretically simulated spectra of GRB 0803 19B integrated over 
the time interval 3s < t£ < 13s with a = — 1.8 (black solid line) and a = 0.0 
(pure thermal case, green dashed line) are compared with the data observed 
by BAT (red points). It can be seen that with the "modified" thermal spectrum 
we can correctly reproduce the observed spectrum, contrary to what happens 
with the pure thermal spectrum. In the residual plot the pure thermal case is 
omitted. 

tra. It must be noted that the upper limit on /p-crb, and there- 
fore the lower limit on B, is less stringent since we cannot 
exclude that more energy has been emitted in the P-GRB out- 
side of the instrumental bandpass. 

Concerning the second episode, lasting from 28 s to the 
end of the prompt emission, we performed numerical sim- 
ulations with different sets of parameters, but we encoun- 
tered several difficulties. In particular, while we can obtain 
a theoretical spectrum compatible with the observed one, it is 
not possible to correctly reproduce the time variability of the 
light curve, even when a bi-dimensional model for the CBM 
is adopted dBianco et al.ll2006t iBernardini etal] |2b09). This 
is consistent with the results presented by other authors: the 
time-resolved prompt emission spectra are best fit with power- 
laws and no change in the photon in dex is observed betwee n 
the first and the second component (Stamat ikos et al.l 12009): 
on the contrary, a va riation of the time-variability is found 
(Marg utti et al.| [2008). A possible explanation for this prob- 
lem is that a fully three-dimensional modeling of the CBM is 
needed. 

5. GRB 050904 

We analyzed the prompt emission light curve (Fig. [9} and 
spectrum (Fig.Qj) of GRB 050904 observed by BAT. As we 
already mentioned in the Introduction, for GRB 050904 we 
have E° b e s ak = 314!'™ keV dSugita et al.H2009h . The data have 
been obtained by using the standard FTOOLS package (Hea- 
soft, version 6.10); the BAT spectru m integrated over the Tgp 
of the source (Tgo = 225 ± 10 s, see ISakamoto et af]|2005l) is 
best modeled with a power-law with photon index 7=1.25 ± 
0.07, with a chi square value of \ 2 = 64.09 for 60 degrees of 
freedom. Within the fireshell model, we identify the prompt 
emission with the peak of the extended afterglow. In this 
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Figure 9. Theoretically simulated light curve of GRB 050904 prompt emis- 
sion in the 15-150 keV energy band (black solid curve) is compared with data 
observed by BAT (red points). The labels "a", "b","c","d" and "e" identify 
the peaks (see Fig.fl0]and Tab.fJ). 
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Figure 10. Structure of the CBM adopted for GRB 050904. The labels "a", 
"b","c","d" and "e" indicate the values corresponding to the peaks in the BAT 
light curve (see Fig.|9]and Tab.|2j. 
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Table 2 

Properties of the CBM structure adopted for GRB 050904: distance from the 

center of the explosion, thickness, normalized density and mass of the 
clumps; for each distance the transverse dimension of the visible area is also 
reported. 
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Figure 11. Theoretically simulated time integrated spectra of GRB 050904 
for < < 225s with a = —1.8 (black solid line) and a = 0.0 (pure thermal 
case, green dashed line) are compared with the data observed by BAT (red 
points). It can be seen that with the "modified" thermal spectrum we can 
correctly reproduce the observed spectrum, contrary to what happens with the 
pure thermal spectrum. In the residual plot the pure thermal case is omitted. 
The ra nge of vertical axes in the residual plot has been chosen to be the same 
of Fig. [ft] 



case the P-GRB has not been observed. In fact, we have 
estimated £S?„ DD = 1.99 x 10 52 erg, that for z = 6.29 corre- 



P-CRB 

sponds to a fiuence of ~ 6.3 x 10~ 9 erg cm -2 . If we assume 
an observed duration Atp-cRB > 1 s, the P-GRB flux is un- 
der the BAT threshold. The numerical simulation that best 
reproduce the observational data is obtained with similar val- 
ues of E\% and B found for GRB 0803 19B: E^ t = 1.0 x 10 54 
erg and B = 2.2 x 10~ 3 , with a Lorentz gamma factor at the 
transparency point 70 = 446. This could be an indication 
of a similar progenitor for the two sources. Concerning the 
other model parameters, we found an average number density 
(n c bm) ~ 0.2 particles cm" 3 and 1Z = 2 x 10~ n ; these values 
are different from the ones obtained for GRB 0803 19B and 
this could be an indication of the fact that the two bursts oc- 
curred in different environments. The structure of the CBM 
adopted is shown in Fig. [10] and the adopted density contrast 
with respect to the mean density is reported in Tab. [2] 

Also in this case the numerical simulation that best repro- 
duces the observational data has been obtained assuming the 
value —1.8 for the free parameter a; in this way we can also 
correctly reproduce spectra integrated over intervals of time 
much less than the of the source (in Fig. [I2]is shown, as 
an example, the BAT spectrum integrated over the first 50 s). 
Once again, the fact that the model is able to reproduce the ob- 
served spectrum regardless of the time scale over which it is 
integrated is a clear support of the correctness of the assumed 
co-moving spectral shape. 

Estimating the sensitivity of the determination of these pa- 
rameters of the model in this case of GRB 050904 is more 
difficult than in the previous case of GRB 0803 19B. In fact, 
in this case the P-GRB is not observed and therefore we 
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Figure 12. Theoretically simulated time integrated spectra of GRB 050904 
for < t$ < 50s with a = — 1.8 (black solid line) and a = 0.0 (pure thermal 
case, green dashed line) are compared with the data observed by BAT (red 
points). It can be seen that with the "modified" thermal spectrum we can 
correctly reproduce the observed spectrum, contrary to what happens with the 
pure thermal spectrum. In the residual plot the pure thermal case is omitted. 

may have only a lower limit on the value of B. However, 
we can proceed as follows. We fix Ef ot to the observed 
value of E^ = 1.04^ x 10 54 erg dSugita et all 120091) . i.e. 
8.7 x 10 53 < Ef* < 1.29 x 10 54 erg. We can then make an 
educated guess about the average value of « c fem> i- e - we can 
assume that 0.1 < (n c b,„) < 10 particles/cm 3 (see Sec. 13.41 ). 
With this choice, we obtain 1.9 x 10" 3 <B<3.4x 10" 3 and 
1.5 x 10" 11 <TZ<8.0x 10" 11 to reproduce the observed light 
curves and spectra. 

6. DISCUSSION ON THE COMOVING SPECTRUM 

We have mentioned in Sec.[T]that GRBs with £■,,„ up to 10 53 
erg have been successfully interpreted within the traditional 
fireshell model by assuming a pure comoving thermal spec- 
trum. In the previous sections we have also shown that a mod- 
ification of the comoving spectrum (see Eq. |7]i is needed to 
correctly reproduce the observational data of two of the most 
energetic GRBs, GRB 0803 19B and GRB 050904, within the 
fireshell model: the difficulty in interpreting the BAT data by 
assuming a pure comoving thermal spectrum for these sources 
has been clearly shown in Fig. [8] 

The reasons of this result can be summarized as follows (see 
also Sec.Q}: 

1. The modification of the spectral energy distribution 
given in Eq. [7] does not affect the spectrum near the 
peak, but it affects the low energy Rayleigh-Jeans tail 
of the distribution (see Figs. Q] and 

2. For the sources with £ !SO < 10 53 erg previously anal- 
ysed within the fireshell model, only the spectral region 
around the peak contributes to the emission in the in- 
strument bandpasses. We give two explicit examples in 
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GRB 


031203 (a) 


060607A (b) 


0803 19B (c) 


050904 (d) 


z 


0.106 


3.082 


0.937 


6.29 


Eis„ (erg) 


10 s0 


1.1 x 10 53 


1.34 x 10 54 


1.04 x 10 54 


Keak (keV) 


158 ±51 


478!™ 


1261 ± 65 


2291«f 


Ep* (keV) 


144 ±46 


117$ 


675 ± 22 


314!™ 


Satellite 


INTEGRAL 


Swift 


Swift 


Swift 


Instrument 


IBIS/ISGRI 


BAT 


BAT 


BAT 


AE (keV) 


17-500 


15-150 


15-150 


15-150 



Table 3 

Redshift (z), isotropic energy (Ej so ), intrinsic and observed peak energy 
(E' peak and E° b e s ak respectively) of GRB spectra analysed within the fireshell 
model. The instruments by which the observational data interpreted within 

the fireshell model have been taken, together with the energy range they 
co ver (AE), are also reported. References for z: (a ) Prochaska et al.ll !2004l) ; 

(b) ILedoux et~aTl <2006l); (c)IVreeswiik et aT] ffflM : Tdl lKawai et aljl2005l) . 
References for £, so : (a) Amati 1 2006); (b) L. Amati, private communication; 

(c) IGolenetskii et aljf2008h ; (d) ISugita etaLl <2009l) . Refe rences for E' peak : 
(a) IAmatil 120061) ; (b) L. Amati, private communication; (c) Golenetskii et al. 

<2TjOH) ; (d) ISugitaet an<2TxBI) . References for E°^ k : (a)IUlanov et al] 
120051) ; (b) L. Amati, priv ate communication; (c) IRacusin et aljj20O8l) ; (d) 
ISugita et afHSTO . 

Tab. [3] the INTEGRAL and Swift BAT observations of 
GRB 031203 and GRB 060607 A only cover the region 
around the peak. 

3. In the case of the most energetic sources, like GRB 
0803 19B and GRB 050904, Ef/ ak is outside the Swift 
BAT bandpass (see Tab. [3]). Conversely, the Swift BAT 
data cover precisely the low energy component of the 
spectrum, where the effect of the new parametrization 
of the spectrum is maximised (see Tab.|3]l. 

The issue if the modification of the comoving spec- 
trum given in Eq. [7J is really universal and applies as well 
to lower energetic sources is still open. It could only 
be settled by future space missions dedicated to the ob- 
servati on of the prompt emission below 10 keV such as 
LOFT dFeroci & The LOFT Consortium! 1201 11) and MIR AX 
d Amati etal .1120 111) . 

As correctly pointed out to us by an anonymous referee, 
a critical test for the modified co moving thermal s pectrum 
can also come from GRB 061121 (Pag e et alj|2007l) . In this 
special case, a soft precursor pulse triggere d BAT, which al- 
lowed Swift to slew in time for BAT, XRT dPage et al.ll2"006h 
and UVOT to simultaneou sly observe the promp t emission. 
Furthermore, Konus-Wind (Golenetskii et al. 200fj) also ob- 
served this burst, which resulted in an E° b e s ak ~ 606 + 90/ -72 

keV (well above the BAT band pass) and Ej SO ~ 2 .5 x 10 53 erg, 
given a redshift of z = 1.314 (Bloom et alj|2006). This source 
presents many other interesting and challenging observational 
features, and a complete analysis is going to be presented in a 
separate paper. However, we like to emphasize that the pho- 
ton index of the XRT data during the main event of the prompt 
emission is 0.6 < T < 0.8 fl in very good agreement with the 
expected value from the modified thermal spectrum assuming 
a = -1 .8, which is indeed ~ 0.8. 

7. CONCLUSIONS 

GRB 050904 was discovered in the pre-Fermi and pre- 
AGILE era, while GRB 0803 19B was discovered in the pre- 

7 see e.g. data at http : / /www . swift . ac . uk/burst analyser/0 



Fermi era and was unobservable by AGILE due to Earth oc- 
cultation. With the exception of the data from Konus- WIND 
and Suzaku-WAM, no observations on the high energy com- 
ponent of these two sources are available. The high quality 
data from the BAT instrument on board the Swift satellite, nev- 
ertheless, have allowed us to reach a quite firm conclusion on 
the low energy component of the spectra of these sources. 

Thanks to these most energetic sources we have been 
able, for the first time, to explore the Rayleigh-Jeans tail 
of the comoving blackbody spectral energy distribution and 
to conclude that it must be modified with an additional 
component. We recall, in fact, that in the original pro- 
posal the thermal nature of the spectrum in the comoving 
frame was adopted only for simplicity, inspired by a simi- 
lar approach followed by Enrico Fermi in the different con- 
text of ultra high-energy collisions. Notice that even pho- 
tospheric emission in ultrarelativisticall y expanding sources 
does not produce pure therm al spectra dRuffini et al.11201 fbl 
see also iPe'er & Rvdell201 ll). The successful interpretation 
of many sources (see e.g. l Bernardini et al. 2005;Ruffini etal. 
2006tlDainotti et alj|2007t iBernardini et al.ll2007t ICaito et all 
2009ll2010Ude Barros et al.1201 lb showed the viability of this 



ansatz (Ru ffini et alj|2004l) . In the intervening years, thanks 
to the data analysis by Felix R yde and collaborators (see e.g. 
lRvddl2004tlRvde & Ps'ei\200% , it has become clear that the 
existence of a pure black body spectra in any G RB o bserva- 
tion is more an exception than the rule (see Sec. 12. H . In the 
present work we show that this is also the case for the co- 
moving spectrum of the extended afterglow. The most inter- 
esting aspect is that it is possible to generalize the previous 
ansatz by the addition of a single power-law component in 
order to recover a consistent interpretation of all previous re- 
sults and of the ones corresponding to the present more en- 
ergetic sources. We have introduced a new phenomenolog- 
ical parameter a describing such an additional component. 
The choice of a = -1.8 leads to a coherent description of 
both sources, not contradictory with the previous results on 
the less energetic sources. The main goal of our work in 
this paper is to maximize the knowledge acquirable for these 
very energetic sources (Ej S0 ~ 10 54 erg), which have a peak 
of emission at energies much higher than the less energetic 
ones (E Pt i > 1 MeV), and therefore to explore the low energy 
part of the spectrum of the prompt emission. A strong and 
promising theoretical activity is currently devoted to ascertain 
a possible role of collisionless shocks in generating power law 
components in the high energy part of the photon spectrum 
(Spitkovskv 2008b). Also, the synchrotr on "line of death" ap- 
pears to be problematic for such models (Sironi & Spitkovskv 
2009). The knowledge of the spectrum in the comoving frame 
is certainly an important step toward the identification of the 
physical process occurring in the interaction of the accelerated 
baryons with the CBM, which is yet largely unknown. 

As an additional result, the analysis of GRB light curves 
and spectra within the fireshell model allows us to infer 
the filamentary, clumpy and porosity structure of the CBM. 
Specifically, we determined (n c b m ) = 6 particles cm" 3 and 
^ = 3.5 x 10- 1() forGRB 080319B, (n cbm )= 0.2 particles cm" 3 
and 11 = 2.0 x 10" 11 for GRB 050904 (see also Tab.©. 

We can correctly reproduce the whole BAT prompt emis- 
sion data of GRB 050904. For GRB 0803 19B only the first 
~ 28 s have been satisfactorily interpreted within the fireshell 
model with a mono-dimensional CBM description. Concern- 
ed 3 9^59 remaining part of the prompt emission, it occurs at 



14 



Patricelli et al. 



GRB 


GRB 0803 19B 


GRB 050904 


E wt (^g) 


1.32 X 10 54 


1.0 x 10 54 


B 


2.3 X 10~ 3 


2.2 x ur 3 


a 


-1.8 


-1.8 


(n c bm) (#/cm 3 ) 


~ 6 


-0.2 


K 


3.5 x 10" 10 


2x 10-" 


Table 4 



Summary of parameters characterizing GRB 0803 19B and GRB 050904. 

a distance r > 10 17 cm, where the transverse dimension of 
the visible area is much larger than the typical size of the 
CBM clumps (5.0 x 10 15 < Ar < 1.7 x 10 16 cm, see Tab.[B. 
Therefore, a fully three-dimensional modeling of the CBM is 
needed at such a distance. 

We have mentioned in Sec. [2] the central role of the phe- 
nomenological parameterizations in the description of GRBs. 
In addition to the a p arameter here defined, we have re called 
the iBand et alj d 19931) formula, the lAmati et all (120021) rela- 
tion and the coefficie nts described by iRydd (12004 120051) and 
iRvde & Pe'erl (|2009). Although no physical explanation for 
these parameters have been reached, they represent certainly 
a fundamental step in reaching a quantitative and qualitative 
description of the source and help to understand the underly- 
ing physical process of GRBs. 

In the fireshell model the radiation observed in the BAT data 
comes from an integration which takes into account the CBM 
filamentary structure and applies a double convolution, over 
the EQTS and the observation time, of a mixing of the co- 
moving thermal and power-law components given by Eq.(|7}. 
In the BAT energy range the cutoffs of the thermal compo- 
nents give a fundamental contribution. It was therefore unex- 
pected that the theoretically computed spectrum would have 
given rise to a power-law so closely resembling the observed 
one. There is no simple relation between the power-law in- 
dex of the observed BAT spectrum and the one of the power- 
law component in the co-moving spectrum given in Eq.©. 
When the co-moving thermal component becomes negligible 
at low enough energy, the co-moving spectrum is described 
by just the power-law component. The convolution of power- 
laws with the same index results in a power-law of that in- 
dex. Observations by XRT in the prompt emission of a highly 
energetic source, if available, may then give direct and in- 
dependent information about the existence of the power-law 
component in Eq.(0 and on its index. If so confirmed, this 
power-law component would not be just a mere phenomeno- 
logical optimization of the agreement between our theory and 
the observed BAT spectra. It would be an independent phys- 
ical component of the co-moving spectrum, whose index can 
be directly read from the observational data. This would give 
an additional strong confirmation of our model. 

We thank C. Guidorzi for the reduced Swift BAT data of 
GRB 0803 19B. We are especially grateful to an anonymous 
referee for her/his important remarks which have improved 
the presentation of our results. 

REFERENCES 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 706, L138 
Ackermann, M., Asano, K., Atwood, W. B., et al. 2010, ApJ, 716, 1178 
Ackermann, M., Ajello, M., Asano, K., et al. 2011, ApJ, 729, 114 
Aksenov, A., Ruffini, R., & Vereshchagin, G. 2007, Phys. Rev. Lett., 99, 
125003 



Aksenov, A. G., Bemardini, M. G., Bianco, C. L., et al. 2010, in SIF 
Conference Proceedings, Vol. 102, The Shocking Universe, ed. 
G. Chincarini, P. D'Avanzo, R. Margutti, & R. Salvaterra, 451 

Amati, L. 2006, MNRAS, 372, 233 

Amati, L., Frontera, F., & Guidorzi, C. 2009, A&A, 508, 173 

Amati, L., Guidorzi, C, Frontera, F, et al. 2008, MNRAS, 391, 577 

Amati, L., Frontera, E, Tavani, M., et al. 2002, A&A, 390, 81 

Amati, L., Feroci, M., Frontera, E, et al. 2011, Nuovo Cimento C, 34, 49 

Aptekar, R. L., Frederiks, D. D., Golenetskii, S. V., et al. 1995, Space 

Science Reviews, 71, 265 
Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 
Band, D., Matteson, J., Ford, L., et al. 1993, ApJ, 413, 281 
Bartolini, C, Greco, G., Guamieri, A., et al. 2009, arXiv:0906.4144 
Bernardini, M. G, Bianco, C. L., Caito, L., et al. 2007, A&A, 474, L13 
Bernardini, M. G., Bianco, C. L., Caito, L., et al. 2008, in American Institute 
of Physics Conference Series, Vol. 1065, 2008 Nanjing Gamma-Ray 
Burst Conference, ed. Y.-F. Huang, Z.-G. Dai, & B. Zhang, 227-230 
Bernardini, M. G, Bianco, C. L., Chardonnet, P., et al. 2005, ApJ, 634, L29 
Bernardini, M. G, Dainotti, M. G., Bianco, C. L., et al. 2009, in American 
Institute of Physics Conference Series, Vol. 1111, Probing Stellar 
Populations out to the Distant Universe, ed. G. Giobbi, A. Tornambe, 
G. Raimondo, M. Limongi, L. A. Antonelli, N. Menci, & E. Brocato, 
383-386 

Bianco, C. L., Caito, L., & Ruffini, R. 2006, Nuovo Cimento B, 121, 1441 
Bianco, C. L., & Ruffini, R. 2004, ApJ, 605, LI 
— . 2005a, ApJ, 633, L13 
— . 2005b, ApJ, 620, L23 

Blandford, R. D., & McKee, C. F. 1976, Physics of Fluids, 19, 1130 
Bloom, J. S., Perley, D. A., & Chen, H. W. 2006, GCN Circ, 5826 
Boer, M., Atteia, J. L., Damerdji, Y., et al. 2006, ApJ, 638, L71 
Caito, L., Amati, L., Bernardini, M. G., et al. 2010, A&A, 521, A80 
Caito, L., Bemardini, M. G, Bianco, C. L., et al. 2009, A&A, 498, 501 
Cavallo, G., & Rees, M. J. 1978, MNRAS, 183, 359 
Crider, A., Liang, E. P., & Preece, R. D. 1998, in American Institute of 
Physics Conference Series, Vol. 428, Gamma-Ray Bursts, 4th Hunstville 
Symposium, ed. C. A. Meegan, R. D. Preece, & T. M. Koshut, 359-363 
Crider, A., Liang, E. P., Smith, I. A., et al. 1997, ApJ, 479, L39 
Cummings, J., Angelini, L., Barthelmy, S., et al. 2005, GCN Circ, 3910, 1 
Curran, P. A., van der Horst, A. J., & Wijers, R. A. M. J. 2008, MNRAS, 
386, 859 

Cusumano, G., Mangano, V, Chincarini, G et al. 2007, A&A, 462, 73 
Daigne, E, Bosnjak, Z., & Dubus, G. 2009, arXiv:0912.3743 
Daigne, E, & Mochkovitch, R. 2002, MNRAS, 336, 1271 
Dainotti, M. G., Bernardini, M. G., Bianco, C. L., et al. 2007, A&A, 471, 
L29 

Damour, T, & Ruffini, R. 1975, Physical Review Letters, 35, 463 

de Barros, G., Amati, L., Bernardini, M. G, et al. 201 1, A&A, 529, A130 

Dermer, C. D. 2006, Nuovo Cimento B, 121, 1331 

— . 2008, ApJ, 684, 430 

Dermer, C. D., & Mitman, K. E. 1999, ApJ, 513, L5 
Eichler, D., & Levinson, A. 2000, ApJ, 529, 146 
Fermi, E. 1949, Physical Review, 75, 1169 

— . 1954, ApJ, 119, 1 

Feroci, M., & The LOFT Consortium. 201 1, arXiv: 1 107.0436 
Finkelstein, S. L., Rhoads, J. E., Malhotra, S., & Grogin, N. 2009, ApJ, 691, 
465 

Frederiks, D., & Pal'Shin, V. 2011, GCN Circ, 12370, 1 
Frontera, F, Amati, L., Costa, E., et al. 2000, ApJSS, 127, 59 
Gehrels, N, Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 
Gendre, B., Galli, A., Corsi, A., et al. 2007, A&A, 462, 565 
Ghirlanda, G., Celotti, A., & Ghisellini, G. 2002, A&A, 393, 409 
— . 2003, A&A, 406, 879 

Ghisellini, G., & Celotti, A. 1999, A&AS, 138, 527 
Giannios, D. 2006, A&A, 457, 763 

Golenetskii, S., Aptekar, R., Mazets, E., et al. 2006, GCN Circ, 5837 

— . 2008, GCN Circ, 7482 

Goodman, J. 1986, ApJ, 308, L47 

Granot, J., Piran, T, & Sail R. 1999, ApJ, 513, 679 

Grazinov, A., & Waxman, E. 1999, ApJ, 511, 852 

Guidorzi, C, Lacapra, M., Frontera, F, et al. 2011, A&A, 526, A49 

Izzo, L., Bernardini, M. G., Bianco, C. L., et al. 2010, Journal of Korean 

Physical Society, 57, 551 

Izzo, L., Ruffini, R., Penacchioni, A. V, et al. 2012. larXiv: 1202.43741 
Kaneko, Y, Gonzalez, M. M., Preece, R. D., Dingus, B. L., & Briggs, M. S. 

2008, ApJ, 677, 1168 
Kaneko, Y, Preece, R. D., Briggs, M. S., et al. 2006, ApJSS, 166, 298 
Kann, D. A., Masetti, N., & Klose, S. 2007, AJ, 133, 1187 



Analysis of GRB 0803 19B and GRB 050904 within the fireshell model: evidence for a broader SED 



15 



Kann, D. A., Klose, S., Zhang, B., et al. 2010, ApJ, 720, 1513 
Kawai, N., Yamada, T., Kosugi, G., Hattori, T., & Aoki, K. 2005, GCN 
Circ, 3937 

Kim, S., Staveley-Smith, L., Sault, R. J., et al. 1998, Pub. Astron. Soc. 

Austr., 15, 132 
Kumar, P. 1999, ApJ, 523, L113 
Kumar, P., & McMahon, E. 2008, MNRAS, 384, 33 
Kumar, P., & Narayan, R. 2009, MNRAS, 395, 472 
Kumar, P., & Panaitescu, A. 2008, MNRAS, 391, L19 
Lazzati, D., & Begelman, M. C. 2010, ApJ, 725, 1137 
Lazzati, D., Ghisellini, G, Celotti, A., & Rees, M. J. 2000, ApJ, 529, L17 
Ledoux, C, Vreeswijk, P., Smette, A., Jaunsen, A., & Kaufer, A. 2006, 

GCN Circ, 5237 
Lockman, F. J. 2002, ApJ, 580, L47 

Margutti, R., Bernardini, M. G, Barniol Duran, R., et al. 201 la, MNRAS, 
410, 1064 

Margutti, R., Guidorzi, C, Chincarini, G, et al. 2008, in American Institute 
of Physics Conference Series, Vol. 1065, 2008 Nanjing Gamma-Ray 
Burst Conference, ed. Y.-F. Huang, Z.-G. Dai, & B. Zhang, 259-262 

Margutti, R., Chincarini, G, Granot, J., et al. 2011b, MNRAS, 417, 2144 

Medvedev, M. V. 2000, ApJ, 540, 704 

Medvedev, M. V, & Loeb, A. 1999, ApJ, 526, 697 

Medvedev, M. V., & Spitkovsky, A. 2009, ApJ, 700, 956 

Meegan, C, Lichti, G, Bhat, P. N, et al. 2009, ApJ, 702, 791 

Meszaros, P. 2002, ARAA, 40, 137 

Meszaros, P., Laguna, P., & Rees, M. J. 1993, ApJ, 415, 181 
Meszaros, P., Ramirez-Ruiz, E., Rees, M. J., & Zhang, B. 2002, ApJ, 578, 
812 

Meszaros, P., & Rees, M. J. 2000, ApJ, 530, 292 

Molinari, E., Vergani, S. D., Malesani, D., et al. 2007, A&A, 469, L13 

Nakar, E., & Granot, J. 2007, MNRAS, 380, 1744 

Paczynski, B. 1986, ApJ, 308, L43 

Page, K. L., Starling, R. L. C, Osborne, J. P., Troja, E., & Morris, D. 2006, 
GCN Circ, 5832 

Page, K. L., Willingale, R., Osborne, J. P., et al. 2007, ApJ, 663, 1 125 
Panaitescu, A., & Meszaros, P. 1998, ApJ, 493, L31 
Panaitescu, A., & Meszaros, P. 2000, ApJ, 544, L17 
Pe'er, A., & Ryde, F. 201 1, ApJ, 732, 49 
Pe'er, A., & Zhang, B. 2006, ApJ, 653, 454 

Penacchioni, A. V., Ruffini, R., Izzo, L., et al. 2012, A&A, 538, A58 

Piran, T. 1999, Phys. Rep., 314, 575 

— . 2004, Reviews of Modern Physics, 76, 1143 

Piran, T., Sari, R., & Zou, Y. 2009, MNRAS, 393, 1 107 

Preece, R. D., Briggs, M. S., Giblin, T. W., et al. 2002, ApJ, 581, 1248 

Preece, R. D., Briggs, M. S., Mallozzi, R. S., et al. 2000, ApJSS, 126, 19 

Preece, R. D., Pendleton, G. N., Briggs, M. S., et al. 1998, ApJ, 496, 849 

Prochaska, J. X., Bloom, J. S., Chen, H.-W., et al. 2004, ApJ, 61 1, 200 

Racusin, J. L., Karpov, S. V., Sokolowski, M., et al. 2008, Nature, 455, 183 

Ramirez-Ruiz, E., & Fenimore, E. E. 2000, ApJ, 539, 712 

Rees, M. J., & Meszaros, P. 1992, MNRAS, 258, 41P 

— . 1994, ApJ, 430, L93 

— . 1998, ApJ, 496, LI 

Ruffini, R., Bernardini, M. G, Bianco, C. L., et al. 2006, ApJ, 645, L109 
Ruffini, R., Bianco, C. L., Chardonnet, P., et al. 2004, IJMPD, 13, 843 
— . 2005, IJMPD, 14, 97 

Ruffini, R., Bianco, C. L., Chardonnet, P., et al. 2003, in American Institute 
of Physics Conference Series, Vol. 668, Cosmology and Gravitation, ed. 
M. Novello & S. E. Perez Bergliaffa, 16-107 

Ruffini, R., Bianco, C. L., Chardonnet, P., Fraschetti, F, & Xue, S.-S. 2001a, 
ApJ, 555, LI 13 



— . 2002, ApJ, 581, L19 

Ruffini, R., Bianco, C. L., Fraschetti, P. C. F, & Xue, S.-S. 2001b, Nuovo 

Cimento B Serie, 1 16, 99 
Ruffini, R., Chakrabarti, S. K., & Izzo, L. 2010a, Submitted to Adv. Sp. Res. 
Ruffini, R., Izzo, L., Penacchioni, A. V, et al. 2011a, PoS(Texas2010), 101 
Ruffini, R., Salmonson, J. D., Wilson, J. R., & Xue, S.-S. 1999, A&A, 350, 

334 

— . 2000, A&A, 359, 855 

Ruffini, R., Siutsou, I. A., & Vereshchagin, G. V. 201 lb. larXivilTT 0.0407 
Ruffini, R., Vereshchagin, G., & Xue, S. 2010b, Phys. Rep., 487, 1 
Ruffini, R., Aksenov, A. G, Bernardini, M. G, et al. 2008, in American 
Institute of Physics Conference Series, Vol. 1065, 2008 Nanjing 
Gamma-Ray Burst Conference, ed. Y.-F. Huang, Z.-G. Dai, & B. Zhang, 
219-222 

Ruffini, R., Aksenov, A. G, Bernardini, M. G, et al. 2009, in American 
Institute of Physics Conference Series, Vol. 1132, XIII Brazilian School 
on Cosmology and Gravitation, ed. M. Novello & S. Perez Bergliaffa, 
199-266 

Ryde, F. 2004, ApJ, 614, 827 

— . 2005, ApJ, 625, L95 

Ryde, F, & Pe'er, A. 2009, ApJ, 702, 1211 

Ryde, F, Axelsson, M., Zhang, B. B., et al. 2010, ApJ, 709, L172 

Ryde, F, Pe'Er, A., Nymark, T., et al. 201 1, MNRAS, 935 

Rykoff, E. S., Aharonian, F, Akerlof, C. W., et al. 2009, ApJ, 702, 489 

Sakamoto, T., Barbier, L., Barthelmy, S., et al. 2005, GCN Circ, 3938 

Sari, R. 1997, ApJ, 489, L37 

— . 1998, ApJ, 494, L49 

Sari, R., & Piran, T. 1999, ApJ, 520, 641 

Schaefer, B. E., Palmer, D., Dingus, B. L., et al. 1998, ApJ, 492, 696 

Shemi, A. 1994, MNRAS, 269, 1112 

Sironi, L., & Spitkovsky, A. 2009, ApJ, 707, L92 

Spitkovsky, A. 2008a, ApJ, 673, L39 

— . 2008b, ApJ, 682, L5 

Stamatikos, M., Ukwatta, T. N„ Sakamoto, T., et al. 2009, in American 
Institute of Physics Conference Series, Vol. 1133, GAMMA-RAY 
BURST: Sixth Huntsville Symposium, ed. C. Meegan, C. Kouveliotou, & 
N. Gehrels, 356-361 
Stern, B. E., & Poutanen, J. 2004, MNRAS, 352, L35 
Sugita, S., Yamaoka, K, Ohno, M., et al. 2009, PASJ, 61, 521 
Tagliaferri, G, Antonelli, L. A., Chincarini, G, et al. 2005, A&A, 443, LI 
Tanvir, N. R., Rol, E., Levan, A. J., et al. 2010, ApJ, 725, 625 
Tavani, M. 1996, ApJ, 466, 768 

Tavani, M., Barbiellini, G, Argan, A., et al. 2009, A&A, 502, 995 
Ulanov, M. V, Golenetskii, S. V, Frederiks, D. D., et al. 2005, Nuovo 

Cimento C, 28, 351 
van Paradijs, J., Kouveliotou, C, & Wijers, R. A. M. J. 2000, ARAA, 38, 

379 

Vreeswijk, P. M., Smette, A., Malesani, D., et al. 2008, GCN Circ, 7444 
Waxman, E. 1997, ApJ, 491, L19 

Wozniak, P. R., Vestrand, W. T., Panaitescu, A. D., et al. 2009, ApJ, 691, 495 
Yamaoka, K, Endo, A., Enoto, T., et al. 2009, PASJ, 61, 35 
Zdziarski, A. A., Svensson, R., & Paczynski, B. 1991, ApJ, 366, 343 
Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354 
Zou, Y.-C, Fan, Y.-Z., & Piran, T. 201 1, ApJ, 726, L2 
Zou, Y.-C, & Piran, T. 2010, MNRAS, 402, 1854 



